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Cd ■ Abstract 

We study localization and chirality properties of eigenvectors of 
jjT^- the lattice Dirac operator. In particular we focus on the dependence 

f-H | of our observables on the size of the corresponding eigenvalue, which 

allows us to study the transition of a near-zero mode into a bulk mode. 
We analyze ensembles of quenched SU(3) configurations using a Dirac 
S^ ■ operator which is a systematic expansion in path length of a solution of 

the Ginsparg- Wilson equation. Our results support the interpretation 
of the excitations relevant for chiral symmetry breaking as interacting 
instantons and anti-instantons. 
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1 Introduction 

Chiral symmetry breaking is one of the most intriguing features of QCD and 
the existence of a non-vanishing chiral condensate is a cornerstone in our 
understanding of hadron phenomenology. Analyzing and understanding the 
mechanisms which lead to the formation of the chiral condensate thus can 
provide deep insights into the nature of relevant excitations of QCD. 

In the last 25 years a phenomenological picture of chiral symmetry break- 
ing through the interaction of instantons and anti-instantons has been de- 
veloped (see e.g. fjl| for recent reviews). A single classical instanton or anti- 
instanton leads to a zero mode of the Dirac operator. An infinitely separated 
instanton anti- instanton pair has two zero modes, but as soon as the fermion 
wave functions start to overlap they mix, and the Dirac operator acquires a 
pair of small complex conjugate eigenvalues. As more instantons and anti- 
instantons are added they build up a non-vanishing density of eigenvalues 
near the origin. This non-vanishing density of eigenvalues in turn is related 
to the chiral condensate via the Banks-Casher relation [g]. 

This picture of chiral symmetry breaking through instantons provides 
a wealth of signatures for the properties of the eigenvectors of the Dirac 
operator with eigenvalues close to 0, the so-called near-zero modes. In par- 
ticular they should have pronounced localization and chirality properties. 
These signatures can be directly tested in an ab-initio calculation using lat- 
tice gauge theory. The recent progress in understanding chiral symmetry on 
the lattice based on the Ginsparg- Wilson relation now allows us to probe 
the QCD vacuum with a Dirac operator which respects chiral symmetry. 
Analyzing properties of eigenvectors of the Dirac operator to learn about 
instantons is a more direct approach than analyzing the gauge fields with 
cooling methods. The eigenvectors directly reflect the relevant excitations 
as seen by the Dirac operator without introducing a filtering process by 
hand. Further evidence for the dominance of instantons has been obtained 
by showing that fermionic n-point functions can be saturated using spectral 
decomposition with only the lowest eigenmodes of the lattice Dirac operator 

III 

Let us stress that while the instanton anti-instanton scenario of chiral 
symmetry breaking is highly popular, it is still far from being established. 
For an alternative based on monopole physics see e.g. M] and references 
therein. Thus detailed numerical tests based on chirally improved lattice 
calculations are highly appropriate. 

In this article we present our results for locality and chirality properties of 
eigenvectors of the lattice Dirac operator. In particular we also focus on the 



dependence of our observables on the size of the corresponding eigenvalue. 
This allows us to study how localization and chirality of a near-zero mode 
change as the corresponding eigenvalue moves into the bulk of the spectrum. 
We analyze ensembles of quenched SU(3) configurations for different lattice 
sizes and three values of the gauge coupling. We use a Dirac operator which 
is a systematic approximation of a solution of the Ginsparg- Wilson equation. 
It has very good chiral properties and is numerically cheaper than an exact 
solution of the Ginsparg- Wilson equation such as e.g. the overlap operator. 
Our results clearly support the interpretation of the excitations relevant for 
chiral symmetry breaking as interacting instantons and anti-instantons. 

The article is organized as follows: In the next section we collect tech- 
nical preliminaries, in particular we discuss our gauge ensembles, the Dirac 
operator we use and the details of the diagonalization. This is followed by 
sections where we discuss our results for the spectral density (Section 3) 
of the Dirac operator and for the localization properties of its eigenvectors 
(Section 4). In Section 5 we analyze the chirality properties of the near- zero 
modes. The article closes with a summary of our results. 

2 Technical remarks 

Gauge configurations: We use ensembles of SU(3) gauge fields in the 
quenched approximation generated with the Liischer-Weisz action || || . In 
addition to the standard plaquette term with coefficient (5\ this action also 
contains a rectangle term with coefficient /?2 and a parallelogram term with 
coefficient /?3 (see || for a detailed discussion of the action) . The coefficient 
Pi is an independent parameter, while /% and (3% can be computed from 
j3\ using e.g. tadpole improved perturbation theory || [HJ. We use [i\ = 
8.10, 8.30 and 8.45. The corresponding values of /% and /?3 were determined 
in [11]. We performed runs on lattices of size 8 4 , 12 4 and 16 4 and the 



statistics for our ensembles is given in Table |lj. We use a mix of Metropolis 
and over-relaxation steps to update the gauge fields. In order to set the 
scale we computed the Sommer parameter [12] r$ for the 16 4 ensembles. 
We give our results for r$/a together with the lattice spacing a (assuming 
tq = 0.5 fm) in Table [j]. 

Dirac operator: Our Dirac operator D is a systematic expansion of a 
solution of the Ginsparg- Wilson equation § , 

75 £ + D 75 = D 75 D. (1) 



r /a 

a 



12 4 
16 4 



j3i = 8.10 /3i 



.30 



3.94(16) 
0.127(5) fm 



4.67(13) 
0.107(3) fm 



Pi = 8.45 

5.00(5) 
0.100(1) fm 



800 
400 
200 



800 
400 
200 



800 
400 
200 



Table 1: Sommer Parameter ro, lattice spacing a and statistics for our gauge 
field ensembles. 



It has been understood during the last few years, that solutions of (|l|) can 
be used to implement chiral symmetry on the lattice. The two knownQ 
exact solutions are perfect actions 12, O] and overlap fermions |15|] . In a 
series of recent publications [11, 16, 17, 18] it has been demonstrated that 
for many problems it is sufficient to work with an approximate solution of 
(111) which is considerably cheaper from a numerical point of view. The idea 
[jl8[ is to expand the most general Dirac operator D on the lattice into a 
sum of simple operators. Each of these simple operators is built from a 
set of paths of the same length. The members of these sets are related by 
symmetry transformations as e.g. rotations, reflections etc. Along each path 
one builds a gauge transporter from the link variables and multiplies it with 
a Kronecker symbol with the endpoints of the path as arguments. The most 
general Dirac operator is then obtained by multiplying this object with an 
element of the Clifford algebra and summing all the resulting simple lattice 



operators with some real coefficients [14, 18]. 

Once the expansion of the most general Dirac operator is established 
it can be inserted into the Ginsparg- Wilson equation. The left hand side 
can be evaluated easily by commuting the Clifford algebra elements in the 
expansion with 75. The right hand side amounts to a commutation with 75 
and multiplications of the individual terms in the expansion of D. It can 
be shown [^] that the product of any two terms of our expansion which 
appears on the right hand side of (|l]) gives a term which is already present 
in the full expansion of D, i.e. shows up also on the left hand side of (Q). 
The coefficients in front of the individual terms on both sides of (|l|) have 
to be equal. This results in a system of coupled quadratic equations which 



can be solved numerically. In [16] an approximate solution with 19 terms 
was constructed and analyzed. We use the same solution here. A detailed 



1 We remark that when one adds a fifth dimension also the domain wall fermions provide 
a solution of the Ginsparg- Wilson equation. 



description of the terms as well as the coefficients can be found in the ap- 



pendix of [ 11 1 . Our expansion of a solution of the Ginsparg- Wilson equation 
is an expansion in terms of path length, since the size of the coefficients 
multiplying each term in the expansion decreases as the length of the paths 
increases Jl6[| . 

The symmetries that are implemented and are responsible for grouping 
the paths into the above mentioned sets are lattice translations and lattice 
rotations as well as parity, charge conjugation and 75-hermiticity, where 
75-hermiticity is defined as 

#75 = 75^. (2) 



The property (g) can be seen 14, 18] to correspond to the properties of D 



which are used to prove CPT in the continuum. Eq. (g) has an important 
implication for the 75-matrix element of the eigenvectors of D. Let ip be 
such an eigenvector, i.e. Dip = \i/j. Then a few lines of algebra show that 



m 



ip 75 4> = , unless A is real . (3) 

This equation is the basis for the identification of topological modes. For 
the continuum Dirac operator and for a D which is an exact solution of 
the Ginsparg- Wilson equation one has exact zero modes. Only zero modes 
have a non- vanishing matrix element with 75. Here we are working with 
an approximate solution of the Ginsparg- Wilson equation and we do not 
have exact zero modes. However, eigenvectors for which the corresponding 
eigenvalue has a non-vanishing imaginary part cannot be candidates for zero 
modes, since their matrix elements with 75 vanish exactly. In our approx- 
imation the topological modes show up as eigenvectors with a small real 
eigenvalue. For the real eigenvalues on the 16 4 lattices we find mean values 
of 0.0124, 0.0122 and 0.0095 for ft = 8.10, ft = 8.30 and ft = 8.45 respec- 
tively. The distributions (see |)20|| ) of the real eigenvalues are very narrow 
and we find that 95% of the real eigenvalues lie between and 0.05. The 
matrix elements ip'^ip typically have a value ~ ±0.8. If one added more 
terms in our approximation then the real modes would be closer to and 
the matrix element with 75 closer to ±1. 

Computation of the eigenvalues: For the computation of the eigenval- 
ues and eigenvectors we used the implicitly restarted Arnoldi method [21]. 



For lattice sizes 8 and 12 we computed 50 eigenvalues and the correspond- 
ing eigenvectors for each of our gauge field configurations and 30 eigenvalues 
and eigenvectors for the ensembles on the 16 4 lattice. The search criterion 



for the eigenvalues was their modulus, i.e. we computed eigenvalues in con- 
centric circles around the origin until 50, respectively 30 were found. The 
boundary conditions for the fermions were periodic in the space directions 
and anti-periodic in time. 

3 Density of eigenvalues 

When describing chiral symmetry breaking based on instanton phenomenol- 
ogy, a major ingredient is the Banks-Casher formula [0] which relates the 
density p(X) of eigenvalues A at the origin to the chiral condensate, 

(^) = -7rp(0) V -1 . (4) 

We use a convention for the density p(X) such that p(X) AA gives the number 
of eigenvalues in the interval AA. The second ingredient for chiral symmetry 
breaking through instantons is the above mentioned interaction of pairs 
of instantons and anti-instantons which leads to a pair of small complex 
conjugate eigenvalues. In the chirally broken phase of QCD the abundance 
of interacting instantons and anti-instantons builds up a finite density of 
eigenvalues near the origin thus leading to a non- vanishing chiral condensate. 
This picture provides a wealth of signatures for the eigenvectors of the Dirac 
operator with small eigenvalues, the near-zero modes. 

Before we discuss our tests of some of these characteristic features of the 
near-zero modes we briefly present our results for the density of eigenvalues. 
For the continuum Dirac operator the spectrum is restricted to the imagi- 
nary axis and the definition of the spectral density is straightforward. For a 
solution of the Ginsparg- Wilson equation ([l]) the spectrum is located on a 
circle of radius 1 with center 1 in the complex plane. For our approximate 
solution D of @) the spectrum fluctuates slightly around this circle (com- 
pare e.g. the spectra shown in [O]). The non- vanishing real parts of the 
eigenvalues of an exact or approximate Ginsparg- Wilson fermion require a 
small modification in the definition of the eigenvalue density p(X). We bin 
the eigenvalues with respect to their imaginary part and count the num- 
ber of eigenvalues in each bin. Eigenvalues with vanishing imaginary part, 
i.e. the zero modes, were left out. After dividing the count in each bin by 
the total number of eigenvalues we obtain the histograms for p(X) used in 
Fig. |H below. Since we are only interested in the density p(A) in a region 
of A where the real parts of the eigenvalues on the Ginsparg- Wilson circle 
are small, binning with respect to the imaginary parts of A gives a good 
approximation of the continuum definition of the spectral density. 
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Figure 1: Density of eigenvalues as a function of IniA. The density is nor- 
malized by the inverse volume of the corresponding lattices. Zero modes 
were left out for the evaluation of the density. 



In Fig. [I] we show our results for the density of eigenvalues. We use 
the spectra for f3\ = 8.10 (top plot), /?i = 8.30 (middle plot) and Pi = 8.45 
(bottom plot) and display the results for all three volumes 8 4 (dashed curve) 
12 4 (dotted curve) and 16 4 (full curve). We remark that for the largest lattice 
(16 4 ) the smallest Matsubara frequency has a value of tt/16 = 0.1963, such 
that many of the eigenvalues are considerably below the smallest eigenvalue 
for the free case. We plot the density p(X) normalized by the volume L 4 as a 
function of Im A where, since the curve is symmetric with respect to reflection 
at the origin, we display only positive ImA. One finds that the data for the 
density essentially fall on universal curves. Only near the origin the density 
is more rounded for the smaller lattices in agreement with universal random 
matrix theory predictions [22|. For the same reason also the curve for the 



largest system (16 4 ) still shows a drop near the origin. Up to this finite-size 
effect, the density remains non-zero down to A = 0, thus building up the 
chiral condensate according to the Banks-Casher formula @. We remark 
that at Pi = 8.45 the smallest lattice (8 4 ) already is in the deconfined phase, 



i.e. there the density develops a gap (compare |llj). It is obvious that for 
Pi = 8.10, i.e. the ensemble with the largest value of the gauge coupling the 
density of eigenvalues and thus the condensate is largest while it becomes 
smaller with increasing Pi. 

We remark that in a recent study |2*| of quenched QCD at finite temper- 
ature using the overlap operator, a small nonvanishing density of eigenvalues 
near the origin was found even in the deconfined phase. In our recent paper 



[11 1 with the ultralocal, chirally improved operator no statistically signif- 
icant signal for a nonvanishing density in the deconfined phase was seen. 
This discrepancy might be related to an observation we made recently |24| |. 
For smooth instantons on the lattice the large extent of the overlap oper- 
ator considerably distorts the zero modes for instantons with radius < 2.5 
in lattice units, i.e. objects with a diameter of < 5.0. We thus suspect that 
using the overlap operator on lattices with time extent 4 (as in p3[) may 
lead to serious finite size effects. 

4 Localization properties 

Let us now come to the analysis of the characteristic features of the near- zero 
modes, i.e. the eigenvectors of the Dirac operator with small eigenvalues. In 
order to develop these characteristic features we first discuss the behavior of 
exact zero modes in the continuum. For an instanton configuration in the 
continuum, the Dirac operator has an exact zero mode if) which is localized 



in space-time around the center of the instanton and this localization has a 
radius proportional to that of the underlying instanton. A gauge invariant 
density which displays this localization is obtained by summing \ip(x,d, c)\ 2 
over the Dirac and color indices d and c at each space-time point x individ- 
ually. A lattice discretization of this scalar density p(x) is given by 

P( x ) = £ H x , c i d Y 4>{x,c,d) , (5) 



where now ip(x,c,d) are the entries of the eigenvectors of our lattice Dirac 
operator D. As expected from the analytic result in the continuum, this den- 
sity shows a clear localization at the position of an isolated smooth instanton 



put on the lattice by hand [24, 25]. We remark that since the eigenvectors 



ip of D are normalized one has 

£p(s) = 1- (6) 

X 

A zero mode of the continuum Dirac operator is also an eigenstate of 75. 
Hence it is interesting to also analyze the gauge invariant pseudoscalar den- 
sity p§(x) defined by 

P5(x) = ^2 ip{x,c,d)* (75 W ip(x,c,d') . (7) 

c,d,d' 

When evaluating the corresponding continuum quantity for a zero mode, 
p$(x) differs from its scalar counterpart only by a possible sign, i.e. the scalar 
density equals the pseudoscalar density for an instanton but the two densities 
have opposite sign for an anti-instanton. If one puts a smooth instanton 
anti-instanton pair on the lattice by hand the density of the corresponding 
near-zero mode shows a dipole behavior. Near the instanton peak of the pair 
p$(x) is positive while it is negative near the center of the anti-instanton. 
When analyzing the densities p(x) and ps(x) for thermalized configurations 
one still finds localized structures and also dipole structures which can be 
interpreted as signatures of isolated instantons or instanton anti-instanton 
pairs (see e.g. [11] for plots of examples). Using cooling techniques to identify 



the instantons independently it has been found that the localization of the 
eigenvectors is concentrated at the same region where the cooling procedure 
finds an instanton (see e.g. ||). 

In order to go beyond analyzing the localization properties of the eigen- 
vectors by merely inspecting plots of the scalar density it is convenient to 



introduce the scalar inverse participation ratio /, which is widely used in 
solid state physics, 

I = VY,P(xf- (8) 

X 

V denotes the volume L 4 . From its definition in Eq. (@) it follows that 
p(x) > for all x. Taking into account the normalization Eq. (||) one finds 
that a maximally localized eigenvector which has support on only a single 
site x' must have p(x) = S xx i. Inserting this into the definition (||) for the 
inverse participation ratio one finds that a maximally localized eigenvector 
has I = V. Conversely, a maximally spread eigenvector has p(x) = 1/V for 
all x. In this case, the inverse participation ratio gives a value of / = 1. 
Thus / has a high value for localized states while I is small for delocalized 
ones. An alternative measure for the localization of the eigenvectors, based 
on their self correlation has been studied in fv. 

We now use the inverse participation ratio to analyze the localization of 
the near-zero modes. According to the instanton picture the small eigenval- 
ues come from interacting instantons and anti-instantons. If an instanton 
and an anti-instanton are infinitely far apart, i.e. they do not interact, the 
Dirac operator still sees them as two independent objects and displays two 
eigenvalues 0. As soon as the two objects interact this degeneracy is lifted 
and the two eigenvalues split into a pair of small complex conjugate eigenval- 
ues. When the instanton and the anti-instanton approach further, the two 
eigenvalues move further up (down) in imaginary direction. At the same 
time the instanton and the anti-instanton start to annihilate. The localiza- 
tion is washed out and the corresponding eigenvalues end up in the bulk of 
the spectrum. Thus one expects that for eigenvalues very close to the origin 
the localization is still large since the two partners remain relatively unper- 
turbed. On the other hand eigenvalues further up (down) the imaginary 
axis come from pairs where the partners are already relatively close to each 
other and should be less localized. 

In Fig. [2] we plot our results for the average (I) of the inverse participation 
ratio as a function of ImA. The histograms were obtained by binning ImA 
and computing the average of I for the eigenvectors with eigenvalues in each 
bin. Topological modes, i.e. modes with small real eigenvalues were left out 
in the computation of the histograms^. Since the distribution is symmetric 
with respect to the origin we display the curves only for ImA > 0. We 



2 We remark that the topological modes have values of / varying between 1.8 and up 
to 100, with the highest values of I coming from defects. The distribution is peaked near 
small values similar to the distribution for the near zero modes shown in Fig. pi 
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Figure 2: The average inverse participation ratio (I) of the near-zero modes 
as a function of ImA for lattice size 16 4 . We remark that in the first bin 
real modes were left out. 



show our results for the 16 4 lattice at all three values of 0i. All three 
curves display their maximum at or near the origin, as expected from the 
picture of interacting instantons and anti-instantons. The same behavior 
was observed for the near-zero modes of the staggered Dirac operator in 
||26|| . Furthermore the eigenvectors for the f3\ = 8.10 ensembles have the 
largest average of / near A = while the eigenvectors for (3\ = 8.45 have 
the smallest values. This poses an interesting question: Is it simply a larger 
number of localized states which leads to a larger condensate for (3\ = 8.10, 
or do in addition also the near zero modes themselves become more localized? 
This question can be analyzed by studying the probability distribution of 
the inverse participation ratio. 

In Fig. U we show the probability distribution P(I) of the inverse par- 
ticipation ratio for the near-zero modes (no zero modes were taken into 
account). More precisely, we show the probability distribution for all eigen- 
vectors with eigenvalues obeying < |Im A| < 0.05, i.e. for those modes for 
which the curves in Fig. ^ show a strong j3\ dependence. The distributions 
P(I) are normalized such that their integral over / equals 1. It is obvious 
from Fig. || that for (3i = 8.10 the distribution is shifted towards larger 
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Figure 3: Distribution of the inverse participation ratio for the near-zero 
modes with < I Im A I < 0.05. 



values of /, while it moves towards smaller values of / as /?i is increased. 
Thus not only the average (I) of the inverse participation ratio is larger for 
smaller /3\, but also the probability of finding a localized mode increases 
with smaller (5\. 

So far we have only concentrated on the localization of the near-zero 
modes. According to the instanton picture of chiral symmetry breaking also 
their local chiral behavior should have an interesting signature. In particu- 
lar p§{x) should have positive sign near instanton peaks and negative sign 
near anti-instanton peaks. In the next section we will discuss an observable 
introduced in |27j which tests the local chirality for each lattice point indi- 
vidually. Before we come to our results for the local chirality let us first look 
at an observable which is a global measure for the chirality of an eigenvector. 

Similarly to the scalar inverse participation ratio / of Eq. (||) we can also 
define a measure 1$ for the localization of the pseudoscalar density p§ (x) (see 

h V V P5 (x) 2 . (9) 



E 

X 



Let us discuss the behavior of I5 for a smooth instanton and a smooth 
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Figure 4: The ratio I5/I of pseudoscalar and scalar inverse participation 
ratio as a function of Im A (in physical units) . The data were computed on 
the 16 4 lattices. 



instanton anti-instanton pair. For a zero mode, corresponding to a single 
smooth instanton or anti-instanton, one has 1 = 1$ since for all x the p{x) 
differ from the ps(x) only by a global sign. For an instanton anti-instanton 
pair the pseudoscalar density p$(x) has to undergo a change of sign, since for 
x near the instanton peak ps(x) is positive while it is negative near the center 
of the anti-instanton. This implies that ps(x) has to go through for some x. 
Since (compare the definitions ©,(0)) for each x we have |ps(x)| < p(x), it 
follows from the definitions of / and I§ that for an instanton anti-instanton 
pair one has ^5 < /. As long as the two constituents are relatively remote 
from each other and much of their identity as an independent instanton, 
respectively anti-instanton remains intact one still finds I5 ~ /. However, 
for pairs where the partners are close to each other, one expects values of 
Is considerably smaller than /. Thus one expects to find a relatively large 
ratio of J5/I for modes with eigenvalues very close to the origin and a drop 
of this ratio as | Im A| increases. 

In Fig. |H we show our results for the average (I5/I) as a function of Im A. 
Again we do not take into account zero modes. One clearly sees the drop 
of (I5/I) as Im A increases. This is in agreement with the above discussed 
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picture of instanton anti-instanton pairs where the two partners lose their 
chiral identity as they approach each other. Note that the curves for the 
three different values of j3\ are very similar. 

5 Local chirality 

To further study the nature of the near-zero modes we analyze the local 
chirality observable proposed by Horvath et al. in J27J. In the original arti- 
cle [27] the computation was done with Wilson fermions and the result was 



interpreted as evidence against the instanton picture. In the meantime the 



local chirality variable was reanalyzed in several articles [28] using different 
Dirac operators or a modification of the observable. These latter collab- 
orations all interpret their results in favor of the instanton picture. Here 
we add further supporting evidence by analyzing for the first time how the 
local chirality of a near-zero mode changes as |ImA| of the corresponding 
eigenvalue increases. This allows us to directly see how an instanton and 
an anti-instanton lose their chirality and the corresponding near-zero mode 
turns into a bulk mode. 

Before we present our numerical results we briefly repeat the definition of 
the local chirality variable of p^ ]. In analogy to the densities p{x) and ps(x) 
of Eqs. (|5|),(|7|) we can define local densities p+(x) and p-{x) with positive 
and negative chirality, 

P±{x) = Yl ^(s,c,d)* {P±)d,d> i>(x,c,J) , (10) 

c,d,d> 

where P± denote the projectors onto positive and negative chirality P± = 
(1 ± 75)/2. It is now interesting to analyze locally for each lattice point 
x the ratio p+(x)/p_(x). For a classical instanton in the continuum, the 
corresponding zero mode ip(x) has positive chirality and the density p+(x) 
is positive for all x while p-{x) vanishes everywhere. Thus the ratio of 
the two always gives oo. For an anti-instanton the roles of p+ and p- are 
exchanged and the ratio is always 0. When now analyzing the eigenvectors 
for an interacting instanton anti-instanton pair this property should still 
hold approximately near the peaks. The ratio p + (x)/p-(x) is expected to 
be large for all x near the instanton peak of the pair and small for all x near 
the anti-instanton peak. In a final step Horvath et al. map the two extreme 
values oo and of the ratio p + {x)/p-{x) onto the two values +1,-1 using 
the arctangent. One ends up with the local chirality variable X{x) defined 
as 
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X(x) = - arctan W , : - 1 . (11) 



If the eigenvectors for the near-zero eigenvalues correspond to instanton 
anti-instanton pairs then the values of X{x) should cluster near +1 and —1 
when one selects lattice points x near the peaks of the scalar density p(x). 
One can choose different values for the percentage of points x used for the 
average and here we will present results for cuts of 1%, 6.25% and 12.5%. 
This means that we average over those 1% (6.25%, 12.5%) of all lattice 
points which support the highest peaks of p(x). The smallest cut-off of 1% 
will give the best signature since only the highest peaks which are not so 
much affected by quantum fluctuations contribute. On the other hand for 
such a small cut-off the results are not expected to be very conclusive, since 
typically instanton models have a packing fraction of instantons considerably 
larger than 1% (see Q). As already announced, we study the dependence 
of the chiral density of the eigenvectors on the imaginary part Im A of the 
corresponding eigenvalues. To do so we bin Im A and sum the contributions 
X{x) for each bin individually. In order to compare the different bins we 
normalize X(x) such that the area under the histogram is 1 for each of the 
A-bins. 

In Fig. H we show our results for the (5\ = 8.10 ensemble on the 16 
lattice and Fig. |6| displays the local chirality for (3\ = 8.45. The top plots 
give the result for the 1% cut on the number of supporting lattice points, 
while the middle and bottom plots are for 6.25% and 12.5% respectively. 
We remark that due to the higher density of eigenvalues A at j3\ = 8.10 our 
data which were computed from the eigenvectors corresponding to the 30 
lowest eigenvalues extend only to |ImA| = 0.1 while we reach |ImA| = 0.15 
at 0i = 8.45. It is obvious, that for all plots the modes with eigenvalues 
closest to have the best signal for local chirality, i.e. the curves in the 
first bin always show two well pronounced maxima in the vicinity of ±1 
and a clear valley separating the two peaks. As |ImA| increases, the signal 
for local chirality weakens, i.e. the valley is less pronounced and the height 
of the peaks lowers. This reflects the washing out of the local chirality 
as instantons and anti-instantons approach, as discussed above. The same 
effect has been observed for the global observable (I5/I) as shown in Fig. ||. 
Thus, as a near-zero mode turns into a bulk mode it loses its local chirality. 

When comparing the results for the different cuts on the percentage of 
the supporting lattice points one finds, as expected, that the strongest signal 
for local chirality is obtained for the 1%-cut where only the highest peaks 



14 




6.25% 




12.5% 




Figure 5: Local chirality for the near-zero modes in the (3\ = 8.10 ensemble 
on the 16 4 lattices. We binned the data with respect to the imaginary part 
|Im A | of the eigenvalues and normalized the histograms to 1 for each A-bin 
individually. We use cuts of 1% (top plot), 6.25% (middle plot) and 12.5% 
(bottom plot) for the percentage of supporting lattice points. 
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6.25% 




12.5% 




Figure 6: Local chirality for the near-zero modes as in Fig. || but now for 
j3i = 8.45. 
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are taken into account. However, also for cuts of 6.25% and 12.5% one still 
finds a clear signal for local chirality. Finally when comparing the results for 
the two different values of j3\ we observe in the (3\ = 8.45 ensemble a better 
signal of local chirality for the modes with eigenvalues very close to 0. This 
can be understood by the stronger suppression of quantum fluctuations at 
larger values of f3\ . 

6 Summary 

In this article we have studied localization and chirality properties of eigen- 
vectors of the lattice Dirac operator with small eigenvalues. According to 
the picture of chiral symmetry breaking through instantons these near-zero 
modes should display characteristic features. In particular, since they are 
supposed to trace interacting instantons and anti- instantons, they are ex- 
pected to be relatively localized and their pseudoscalar density should dis- 
play dipole behavior. We analyze these modes using the scalar and pseu- 
doscalar participation ratio and the local chirality variable of Horvath et 
al. In particular we focus on the dependence of these observables on the 
imaginary parts ImA of the corresponding eigenvalues. One expects that 
both the locality and the local chirality of the near-zero modes decrease 
with increasing | ImA |. Our main findings are as follows: 

• When plotting the average inverse participation ratio (I) as a func- 
tion of ImA, we find the largest values of (/) near the origin and {I) 
decreases as | Im A | increases. 

• A comparison of the probability distribution P(I) of the inverse par- 
ticipation ratio for the near-zero modes shows that for the strongest 
gauge coupling (i.e. f3\ = 8.10), where the chiral condensate is biggest, 
one has the largest probability of finding localized states. 

• A plot of the average ratio (I5/I) as a function of ImA shows that the 
most chiral states are located near the origin and as | Im A | is increased 
this chirality becomes washed out. We observe that the ratio (I5/I) 
falls on a universal curve for the different values of f3\, i.e. is relatively 
stable under quantum fluctuations. 

• The observable of Horvath et al. gives a clear signal for local chirality 
for cuts of 1%, 6.25% and 12.5% on the number of supporting lattice 
points. When | Im A | is increased we observe a decay of the local 
chirality for the corresponding modes. 
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Our results support the picture of the emergence of the chiral condensate 
through the interaction of instantons and anti-instantons. 
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